Dirac-like plasmons in honeycomb lattices of metallic nanoparticles 
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We consider a two-dimensional iioneycomb lattice of metallic nanoparticles, each supporting a localized 
surface plasmon, and study the properties of the collective plasmons resulting from the near field dipolar inter- 
action between the nanoparticles. We analytically investigate the dispersion, the effective Hamiltonian and the 
eigenstates of the collective plasmons for an arbitrary orientation of the individual dipole moments. When the 
polarization points close to the normal to the plane the spectrum presents Dirac cones, similar to those present 
in the electronic band structure of graphene. We show that the corresponding eigenstates of the collective plas- 
mons represent Dirac-like massless bosonic excitations. We further discuss how one can manipulate the Dirac 
points in the Brillouin zone and open a gap in the collective plasmon dispersion by modifying the polarization 
of the localized surface plasmons, paving the way for a fully tunable plasmonic analogue of graphene. 
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Light has been the source of inspiration for scientific think- 
ing for millennia. Ancient Assyrians developed the first lenses 
in order to bend the trajectory of light and control its propa- 
gation. In contrast to the macroscopic scale, the use of light 
to observe microscopic structures poses difficulties due to the 
diffraction limit [1], In an attempt to overcome this limit and 
observe subwavelength structures, plasmonic nanostructures 
have been created [2, 3], like isolated metallic nanoparticles 
[4]. The evanescent field at the surface of the nanoparticle, as- 
sociated to the localized surface plasmon (LSP) resonance [5], 
produces strong optical field enhancement in the subwave- 
length region, allowing one to overcome the diffraction limit 
and achieve resolution at the molecular level [6]. 

While the field of plasmonics mostly focuses on single or 
few structures, the creation of ordered arrays of nanoparticles 
constitutes a bridge to the realm of metamaterials. This inter- 
esting contamination resulted in plasmonic metamaterials that 
exhibit unique properties beyond traditional optics, such as 
negative refractive index [7], perfect lensing [8], the exciting 
perspective of electromagnetic invisibility cloaking [9], and 
"trapped rainbow" slow light exploiting the inherent broad- 
band nature of plasmonics [10]. Indeed, in plasmonic meta- 
materials the interaction between LSPs on individual nanopar- 
ticles generates extended plasmonic modes involving all LSPs 
at once [11, 12]. Understanding the nature and properties of 
these plasmonic modes [referred to as "collective plasmons" 
(CPs) in what follows] is of crucial importance as they are the 
channel guiding electromagnetic radiation with strong lateral 
confinement over macroscopic distances. 

CPs in periodic arrays of metallic nanoparticles are an ac- 
tive area of research in plasmonics because the interaction of 
the LSP resonances can lead to dramatic changes in the overall 
optical response of such structures. For example, it was both 
predicted [13] and observed [14] that the plasmonic response 
of a periodic array of nanoparticles could be significantly nar- 
rowed with respect to the single particle response if the in- 



terparticle separation was of the order of the resonant wave- 
length of the LSP. Further work has shown that these coupled 
resonances are relevant to applications in light emission. For 
example, the direction and polarization of light from adjacent 
fluorescent species can be controlled using such arrays [15]. 
The extended nature of the collective resonances means that 
there is scope for harvesting emission from sources spread 
over large volumes. AiTays of metallic nanoparticles may be 
designed to exhibit stop-gaps and band-edges [16] and are also 
being intensively investigated in the context of Ught harvest- 
ing for photovoltaic devices [17]. 

The dispersion of CPs and their physical nature crucially 
depend on the lattice structure of the metamaterial and on the 
microscopic interaction between LSPs. A lattice which re- 
cently generated remarkable interest in the condensed matter 
community is the honeycomb structure exhibited by graphene, 
a two-dimensional (2D) monolayer of carbon atoms [18]. In 
the case of graphene, the hopping of electrons between neigh- 
boring atoms gives rise to a rich band structure character- 
ized by the presence of fermionic massless Dirac quasipar- 
ticles close to zero energy [19]. The chirality associated 
with pseudo-relativistic Dirac fermions results in several of 
the remarkable properties of graphene, such as a nontrivial 
Berry phase accumulated in parallel transport [20] and the 
suppression of electronic backscattering from smooth scat- 
terers [21], which is ultimately responsible for the very high 
mobility of graphene samples. Undoubtedly, it would be ex- 
citing to harvest the remarkable physical properties of elec- 
trons in graphene in suitably designed plasmonic metamateri- 
als ("plasmonic graphene") by analyzing the Hamiltonian and 
the consequent nature of CP eigenmodes in 2D honeycomb 
lattices of metallic nanoparticles. This is the purpose of the 
present theoretical paper 

We analytically show how the problem of interacting LSPs 
in the honeycomb structure can be mapped to the kinetic 
problem of electrons hopping in graphene, yielding massless 
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FIG. 1. (color online), (a) Honeycomb lattice with lattice constant 



, and lattice vectors ai — a(\/3, 0) and a2 = a( 
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The three 



vectors ei = a(0, -1), 62 = a(^, |) and 63 = a(-^, |) 
connect the A and B inequivalent lattice sites (blue/light gray and 
red/dark gray dots in the figure), (b) First Brillouin zone in reciprocal 
space with primitive vectors bi = |f ("s/S, —1) 



andb2 = f (0,1). 



Dirac-like bosonic CPs in the vicinity of two Dirac points in 
the Brillouin zone. While the conical dispersion of plasmons 
in a honeycomb lattice of nanoparticles has been discussed in 
the past for out-of -plane or purely in-plane polarization [22], 
here we unveil the full Hamiltonian of quantum CPs as well 
as the pseudospin structure of the CP eigenmodes for dipolar 
LSPs with arbitrary orientation. The existence of Dirac points 
is robust for a small in-plane component of the polarization, 
where the system maps to strained graphene [23], while band 
gaps can emerge for increasing in-plane polarization. At en- 
ergies away from the Dirac point, van Hove singularities [24] 
emerge in the CP density of states (DOS), associated with Lif- 
shitz transitions in the topology of equipotential lines [25]. 
Our analysis highlights the physical nature of CP eigenmodes 
as well as the tunability of their band structure and of the cor- 
responding DOS with the polarization of light, which can be 
crucial for enhancing the coupling of light with the plasmonic 
metamaterial at different wavelengths. 

Specifically, we consider an ensemble of identical spher- 
ical metallic nanoparticles of radius r forming a 2D honey- 
comb lattice with lattice constant a embedded in a dielec- 
tric medium with dielectric constant (see Fig. 1). The 
nanoparticles are located at positions Rg, with s — A, B a 
sublattice index which distinguishes the inequivalent lattice 
sites. Each individual nanoparticle supports a LSP resonance 
which can be triggered by an oscillating external electric field 
with wavelength A much larger than r. Under such a con- 
dition, the LSP is a dipolar collective electronic excitation at 
the Mie frequency loq — uj-p/^/l + 2eni, which typically lies 
in the visible or near-infrared part of the spectrum [5]. Here, 
Wp = •\/47r7ioe^/mo is the plasma frequency, with Uc, — e and 
me the electron density, charge, and mass, respectively. The 
LSP corresponding to the electronic center-of-mass excitation 
can be generally considered as a quantum bosonic mode, par- 
ticularly when the size of the nanoparticle is such that quanti- 
zation effects are important [26-29]. The noninteracting part 
of the Hamiltonian describing the independent LSPs on the 



honeycomb lattice sites reads [27, 28] 
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where hsCR) is the displacement field associated with the 
electronic center of mass at position R, Hg (R) its conjugated 
momentum and M = Nciric its mass, with the number of 
valence electrons in each nanoparticle. 

The nature of the coupling between LSPs in different 
nanoparticles depends on their size and distance. Provided 
that the wavelength associated with each LSP is much larger 
than the interparticle distance a and that r < a/3 [12], each 
LSP can be considered as a point dipole with dipole mo- 
ment p = — eiVo/is(R)p which interacts with the neighbor- 
ing ones through dipole-dipole interaction. Moreover, it has 
been numerically shown [22] that a quasistatic approximation 
which only takes into account the near field generated by each 
dipole qualitatively reproduces the results of more sophisti- 
cated simulations in which retardation effects are included. 
Within such a quasistatic approximation, the interaction be- 
tween two dipoles p and p' located at R and R', respectively, 
reads [30] V = [p • p' - 3(p • n)(p' • n)]/e^\R - R'p with 
n = (R — R')/|R — R'|. In what follows, we assume that in 
a CP eigenmode all nanoparticles are initially polarized in the 
same direction p = sm9{sm(px — cos(^y) + cosflz, where 
6 is the angle between p and z, and ip the angle between the 
projection of p in the xy plane and ei [see Fig. 1(a)]. We thus 
write the total Hamiltonian of our system of coupled LSPs as 
H = Hq + Hint, where Hq is defined in Eq. (1) and where 
the dipole-dipole interaction term reads 



H, 
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27r[j — l]/3), and the vec- 
tors Bj connect the A and B sublattices [see Fig. 1(a)]. Notice 
that in Eq. (2), we have only considered the dipole-dipole in- 
teraction between nearest neighbors, as the effect of interac- 
tions beyond nearest neighbors does not qualitatively change 
the plasmonic spectrum [31]. 

The analogy between the plasmonic structure of Fig. 1 
and the electronic properties of graphene becomes transpar- 
ent by introducing the bosonic ladder operators aR|6R — 
(Afwo/2fi)^/^/iA|B(R-)+inA|B(R.)/(2M^a;o)i/2 ^hi^-h sat- 
isfy the commutation relations [aj^,a^,] = [fep^jfe^^^/] = <5r^r' 
and [flp^i^R/] = 0. As we will show in the following, 
the introduction of such operators not only gives access to 
the plasmon dispersion (which can be calculated classically 
as well [31]), but also unveils the Dirac nature of the CP 
quantum states. The harmonic Hamiltonian (1) can be writ- 
ten in terms of the above-mentioned bosonic operators as 

Hn = /iwoEra «Ra"Ra + ^^oErb ^LVb' while Eq. 
(2) transforms into 
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In Eq. (3), n = uJo{r/a)^{l + 2ein)/6em, such that ft < ujq. 
The first term on the right-hand side of Eq. (3) (together with 
its hermitian conjugate) resembles the electronic tight-binding 
Hamiltonian of graphene [19], except for three major differ- 
ences: (i) The Hamiltonian of graphene describes fermionic 
particles (electrons), while we deal here with bosonic excita- 
tions (LSPs). (ii) In graphene, an electron "hops" from one 
lattice site to a neighboring one, i.e., the underlying mecha- 
nism linking the two inequivalent sublattices is purely kinetic. 
In the present case, the mechanism coupling the two sublat- 
tices is purely induced by near-field (dipolar) interactions, 
leading to the creation of an LSP excitation at lattice site Rb 
and the annihilation of another LSP at a nearest neighbor lo- 
cated at Rb + Bj . (iii) In (unstrained) graphene, the hopping 
matrix element between two neighboring atoms is the same 
for all three bonds, while in our case, the three energy scales 
MlCj ij — 1, 2, 3) are in general different and can be tuned 
by the direction of the polarization p of the CP eigenmode 
[cf. Eq. (2)]. For < 9 0o and tt ~ Oq 9 < n with 
6*0 = arcsin \/l /3, the coefficients Cj are all positive for any 
(fi and have different values, resulting in different couplings 
between the bonds, thus mimicking the effect of strain in the 
lattice [23]. For 9o < 9 < tt — 9q, the signs of the coefficients 
Cj depend on (p, and the analogy with strained graphene is 
no longer valid. In the special case where Ci = C2 = C3 
(for = or tt), we expect the CP spectrum to resemble that 
of the electronic band structure in graphene, since the Bloch 
theorem does not depend on the quantum statistics of the par- 
ticles one considers, but only on the structure of the periodic 
lattice [32]. As we will now show, two slight differences ap- 
pear in the CP dispersion as compared to the graphene band 
structure, i.e., the effect of the Hamiltonian Hq is to produce 
a global energy shift (by an amount IIujq), while the "anoma- 
lous" term cx ^RgflRg+e ™ ^'i- (•^) introduces corrections of 
order (il/wo)^ to the spectrum. 

Introducing the bosonic operators in momentum space 
Oq and 6q through the Fourier decompositions aR|6R — 
X]q=(g^ g„) "sxp (iq ■ R)aq|foq, with Af the number 
of unit cells of the honeycomb lattice, the Hamiltonian 
H = Hq + iJint transforms into H = tii^o^qio-qO'q + 
bl^b^) + hnj:^[f^b\{a^ + al^) + H.c] with /q - 

6xp (iq • Gj). The latter Hamiltonian is diagonalized 
by two successive Bogoliubov transformations [33]. First, we 
introduce the two bosonic operators 

in terms of which we obtain H — X]r=±Sq[(^o + 

Thn\f^\)a-[^^a^^ + r^^(a^^aLq^ + H.c.)]. Sec- 
ond, we define two new bosonic modes = 

coshz^qOfq — sinh 1?^ 0;^^^ , with cosht9q = 
2-i/2[(i±f]|/q|/c.„)/(l±2r!|/q|/c.o)i/2 + and 
sinh^± = T2-i/2[(l±r!|/q|/wo)/(l±2r!|/q|/c^o)i/' - 
1]^/^, which diagonalize the Hamiltonian H as 



H = J2r=±J2cif^qPq l^q- The two CP branches 
have dispersions 



-WoJl±2— |/q| (5) 

V Wo 

which reduce to uj^ ~ wq ± ri|/q| to first order infl/ujQ <^ 1, 
for which we have — Q^q ■ Notice that the dispersion (5) 
can also be obtained classically by normal mode analysis [31]. 

The dispersion (5) is shown in Fig. 2 in the case of a polar- 
ization p perpendicular to the plane of the honeycomb lat- 
tice [9 — 0, Fig. 2(a)], in the case of an in-plane polar- 
ization [9 = tt/2, (f ~ 0, Fig. 2(b)], and for the special 
case 9 — arcsin (p — [Fig. 2(c)]. In the first 

case [Fig. 2(a)], we have gapless modes with two inequiva- 
lent Dirac cones centered at the K and K' points located at 
±K = ^^(±1, 0) in the first Brillouin zone [cf. Fig. 1(b)], 
while in the second case, the modes are gapped [Fig. 2(b)]. 
The dispersion shown in Fig. 2(c) corresponds to a polariza- 
tion for which Ci = in Eq. (2), i.e., the bonds linked by 
Gi [cf. Fig. 1(a)] are ineffective and the system is effectively 
translationally invariant along one direction. Hence, the CP 
dispersion in Fig. 2(c) does not depend on qy and presents 
Dirac "lines". 

The analogy between the dispersion shown in Fig. 2(a) 
and the electronic band structure of graphene [19] is striking. 
Close to the two inequivalent Dirac points K and K' [see Fig. 
1(b)], the function /q expands as /q ~ ^{T^qx + ^^qy) with 
q = ±K + (5q(|(5q| < |±K|), such that the dispersion (5) is 
linear and forms a Dirac cone, ~ wq ± w|^q|, with group 
velocity v = 3ila/2. This feature is consistent with numeri- 
cal analysis [22]. In the case of graphene, the purely kinetic 
Hamiltonian is iJ = — i /q^q^q + H.c, with flq and 6q 
fermionic operators, and t the nearest-neighbor hopping en- 
ergy (in the case of unstrained graphene, all coefficients Cj 
are equal to one in the function /q). The consequent diago- 
nal Hamiltonian H = ~ ^ Tt\fc^\a^ is associated with 
the very same eigenstates Oq as in Eq. (4), where the ampli- 
tudes of sublattice A and B are linked by ±/q/|/q|. In the 
vicinity of the Dirac point K (K'), the Hamiltonian in the AB 
basis is expanded as H = hvcr ■ 5q (H = hvcr* ■ Sq) with cr 
the vector of Pauli matrices. This massless Dirac Hamiltonian 
is fulfilled also by CPs near the Dirac points (up to a global 
energy shift of Hloq) and the corresponding states are charac- 
terized by chirality cr ■ 6q ~ ±1. As a consequence, CPs 
will show similar effects to electrons in graphene like a Berry 
phase of tt [20] and the absence of backscattering off inho- 
mogeneities [21]. This could have crucial implications for the 
efficient plasmonic propagation in array-based metamaterials. 

In Fig. 2, the panels (d)-(f) show the density of states 
(DOS) corresponding to the spectrum illustrated in the pan- 
els (a)-(c). It is interesting to notice the tunability of the DOS 
with the direction of the polarization, as well as the emergence 
of van Hove singularities [24]. The latter are associated with 
Lifshitz transitions [25] in the topology of equipotential lines 
that percolate at specific energies. The tunability of van Hove 
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FIG. 2. (color online), (a)-(c) Collective plasmon dispersion relation from Eq. (5) and (d)-(f) corresponding density of states for (a),(d) the 
out-of-plane mode (9 = 0), (b),(e) one in-plane mode (6 — tv/2), and (c),(f) 9 ~ arcsin (g) Polarization angles {6, (p) for which the 

collective plasmon dispersion is gapless (dark blue regions) and gapped (white regions), (h) Collective plasmon dispersion along the K'FK 
direction (qy = 0) for different orientations 6 of the dipoles. In the figure, (p = and Q/ujq = 0.01. 



singularities in the spectrum could be of crucial importance 
to increase the coupling of light of different wavelengths with 
extended CP modes. 

For an arbitrary polarization of the LSPs, we can determine 
if the CP dispersion is gapless by imposing |/q| = in Eq. (5), 
which leads to the condition sC [(C2+C3)2-Cf]/4C2C3 ^ 1 
for having gapless plasmonic modes [31]. In Fig. 2(g), we 
show in dark blue the regions of stability of a massless Dirac 
spectrum in the {9, ip) parameter space for which one has gap- 
less plasmon modes, an example of which is shown in Fig. 
2(a). In Fig. 2(g), the white regions correspond to polariza- 
tions for which the CP dispersion is gapped [as an example, 
see Fig. 2(b)]. Thus, changing the polarization allows one to 
qualitatively change the CP spectrum. This is further illus- 
trated in Fig. 2(h) where we show the CP dispersion (5) along 
the KTK direction [see Fig. 1(b)] for different angles 9 of the 
polarization (in the figure, (p — 0). As one can see from Fig. 
2(h), the two inequivalent Dirac points located at K and K' 
for 6* = drift as one increases 9 and they merge at q = for 
9 ~ arcsin y^2/3, forming parabolic bands, to finally open 
a gap for 9 > arcsin \/2/3 (exemplified by 6* = 7r/2 in the 
figure). 

A limitation on the experimental observability of the CP 
dispersion is plasmonic damping, which tends to blur the res- 
onance frequencies. In order to estimate the feasibility of 
such experiments, we compare the bandwidth of the CP dis- 
persion to the losses in individual nanoparticles. In the lat- 
ter, two main sources of dissipation arise: (i) radiation damp- 
ing with decay rate j^g^^ = 2r^ujQ/3c^ [34] (c is the speed 
of light) which dominates for larger nanoparticle sizes, and 
(ii) Landau damping with decay rate 7l — 3vpg/Ar [26, 28] 



[vp = /i(37r^nc)^/'^/mo is the Fermi velocity and g a con- 
stant of the order of one] which dominates for smaller sizes. 
Hence, there exists an optimal size ropt = {Svpgc^ /Sy^'^ /luo 
for which the total damping 7tot — 7rad + 7l is mini- 
mal. For Ag nanoparticles, we find ropt — 8nm for which 
7tot = 0-leV/h. Note that in Ag, the interband transitions 
are decoupled from the LSP, contrary to the case of Au which 
shows a significantly broader resonance [35]. With an in- 
terparticle distance a — Svopt which maximizes the dipolar 
coupling between nanoparticles [12], we find for Ag that the 
bandwidth (for em = 1, and for the out-of-plane polarization) 
is of the order of Auj — ujq — ojq = wo/9 ~ 0.6cV/h at 
the center of the Brillouin zone. Thus Aoj is sufficiently large 
when compared to 7tot that the plasmon excitation is well de- 
fined and hence clearly measurable. We expect this to be true 
even in the presence of inhomogeneous broadening resulting 
from the size dispersion of the nanoparticles. Moreover, the 
appropriate use of active (gain-enhanced) media [36] might 
increase the observability of the CP dispersion and open the 
way towards a new class of Dirac-point lasing. 

A last comment is in order about the excitation of CPs by 
external photons, whose in-plane momentum must match the 
plasmonic one. In fact, the vicinity of the Dirac points typ- 
ically lies outside the light cone. In order to overcome this 
momentum mismatch and observe the Dirac -like plasmons, 
one might add an extra periodic modulation of the lattice to 
allow grating coupling between the incident light and the de- 
sired collective modes [3]. Another alternative might be to use 
a non-hnear technique to overcome the momentum mismatch 
[37]. 

In conclusion, we demonstrated the strong analogies be- 
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tween the physical properties of electrons in graphene and 
those of collective plasmon modes in a 2D honeycomb lat- 
tice of metallic nanoparticles. While we considered spherical 
nanoparticles, our results also extend in principle to arbitrary 
shapes provided they support dipolar excitations. Whereas 
the electronic states of graphene can be described by mass- 
less Dirac fermions, the CP eigenstates correspond to mass- 
less Dirac-like bosonic excitations. The spectrum of the latter 
can be fully tuned by the polarization of an external light field, 
opening exciting new possibilities for controlling the propaga- 
tion of electromagnetic radiation with subwavelength lateral 
confinement in plasmonic metamaterials. 

We thank Cosimo Gorini, Rodolfo Jalabert, and Dietmar 
Weinmann for stimulating discussions. 
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